NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Technical Memorandum 33-590 

Guidance Strategies and Anaiysis for 
Low -Thrust Navigation 

R. A. Jacobson 


p $3.50 CSCL 17g 


N73-1866i 


Onclas 

GJ/z-l 63447 


JET PROPULSION LABORATORY 
CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 



February 1, 1973 



TECHNICAL REPORT STANDARD TITLE PAGE 


1. Report No. 33-590 


4. Title and Subtitle 


2. Government Accession No. 3. Recipient's Catalog No. 


GUIDANCE STRATBSIES AND ANALYSIS FOR 
LOW-THRUST NAVIGATION 


7, Author (s) 


R. A. Jacobson 


9. Performing Organization Name and Address 

JET PROPULSION LABORATORY 
California Institute of Technology 
4800 Oak Grove Drive 
Pasadena, California 91103 


12. Sponsoring Agency Name and Address 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Washington, D.C. 20546 


15. Supplementary Notes 


5. Report Date I 973 


6 . Performing Organization Code 


8 . Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 

NAS 7-100 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


16. Abstract 

A low-thrust guidance algorithm suitable fo r nppr a ti nm1 nn idiWliii— j— 
lated. A constrained linear feedback control law has been obtained using a 
minimum terminal miss criterion and restricting control corrections to 
constant changes for specified time periods. Both fixed- and variable- 
time -of -arrival guidance were considered. The performance of the guidance 
law was evaluated by applying it to the approach phase of the I 98 O rendez- 
vous mission with the comet Encke. 


17. Key Words (Selected by Author(s)) 

Control and Guidance 
Orbits and Trajectories 
Propulsion, Electric 


18. Distribution Statement 

Unclassified — Unlimited 


19. Security Clossif. (of this report) 20. Security Classif. (of this page) 21 . No. of Pages 22. Price 
Ifeclassified Unclassified 22 
















HOW TO FILL OUT THE TECHNICAL REPORT STANDARD TITLE PAGE 


Make items 1, 4, 5, 9, 12, and 13 agree with the corresponding information on the 
report cover. Use all capital letters for title (item 4). Leave items 2, 6, and 14 
blank. Complete the remaining items as follows: 

3. Recipient's Catalog No. Reserved for use by report recipients. 

7. Author(s). Include corresponding information from the report cover. In 
addition, list the affiliation of an author if it differs from that of the 
performing organization. 

8. Performing Organization Report No. Insert if performing organization 
wishes to assign this number. 

10. Work Unit No. Use the agency-wide code (for example, 923-50-10-06-72), 
which uniquely identifies the work unit under which the work was authorized. 
Non-NASA performing organizations will leave this blank. 

1 1. Insert the number of the contract or grant under which the report was 
prepared. 

15. Supplementary Notes. Enter information not included elsewhere but useful, 
such as: Prepared in cooperation with. . . Translation of (or by). . . Presented 
at conference of. . . To be pubi ished in. . . 

(not to exceed 200 words) factual summary of the 
most significant information contained in the report. If possible, the 
abstract of a classified report should be unclassified. If the report contains 
a significant bibliography or literature survey, mention it here. 

17. Key Words. Insert terms or short phrases selected by the author that identify 
the principal subjects covered in the report, and that are sufficiently 
specific and precise to be used for cataloging. 

18. Distribution Statement. Enter one of the authorized statements used to 
denote releasabil ity to the public or a limitation on dissemination for 
reasons other than security of defense information. Authorized statements 
are "Unclassified— Unlimited, " "U. S. Government and Contractors only, " 

"U. S. Government Agencies only, " and "NASA and NASA Contractors only. " 

19. Security Classification (of report). NOTE: Reports carrying a security 
classification will require additional markings giving security and down- 
grading information as specified by the Security Requirements Checklist 
and the DoD Industrial Security Manual (DoD 5220. 22-M). 

20. Security Classification (of this page). NOTE: Because this page may be 
used in preparing announcements, bibliographies, and data banks, it should 
be unclassified if possible. If a classification is required, indicate sepa- 
rately the classification of the title and the abstract by following these items 
with either "(U)" for unclassified, or "(C)" or "(S)" as applicable for 
classified items. 

21. No. of Pages. Insert the number of pages. 

22. Price. Insert the price set by the Clearinghouse for Federal Scientific and 
Technical Information or the Government Printing Office, if known. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Technical Memorandum 33-590 

Guidance Strategies and Anaiysis for 
Low -Thrust Navigation 

R. A. Jacobson 


JET PROPULSION LABORATORY 
CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 


February 1, 1973 



Prepared Under Confract No. NAS 7- 1 00 
National Aeronautics and Space Administration 









PREFACE 


The work described in this report was performed by the 
Analysis Division of the Jet Propulsion Laboratory. 


Mis 


Preceding pege blenk 


JPL Technical Memorandum 33-590 





Not 


CONTENTS 


I. Introduction 1 

II. The Low- Thrust Guidance Algorithm 2 

III. Statistical Guidance Analysis 8 

IV. Guidance Algorithm Evaluation -- Encke Rendezvous Guidance. . . 9 

V, Concluding Remarks 12 

References 13 

Appendix A. The Reduced J Function and the Final Time Variation ... 19 

Appendix B. The Weighting Matrix W 20 


FIGURES 



1. 

RMS final state errors vs flight time and weighting 
parameters, mode 1 guidance 

15 

2. 

RMS final state errors vs flight time and weighting 
parameters, mode 2 guidance 

15 

3. 

RMS final state errors vs weighting parameters for 30-day 
rendezvous, mode 3 guidance 

16 

4. 

RMS final state errors vs weighting parameters for 30-day 
rendezvous, mode 4 guidance 

16 

5. 

RMS final state errors with reduced control constraints for 
30-day rendezvous, mode 2'guidance 

17 

6. 

RMS final state errors with reduced control constraints for 
30-day rendezvous, mode 4 guidance 

18 






JPL Technical Memorandum 33-590 


V 



ABSTRACT 


A low-thrust guidance algorithm suitable for operational use has been 
formulated. A constrained linear feedback control law has been obtained 
using a minimum terminal miss criterion and restricting control corrections 
to constant changes for specified time periods. Both fixed- and variable- 
time - of- arrival guidance were considered. The performance of the guidance 
law was evaluated by applying it to the approach phase of the 1980 rendezvous 
mission with the comet Encke. 


VI 
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I. INTRODUCTION 


Guiding a low-thrust spacecraft requires techniques differing from 
those which have been employed for ballistic spacecraft. The low-thrust 
engine, which provides a continuous, controllable accelerating force acting 
on the vehicle, not only increases its targeting ability but also increases the 
demands imposed upon its navigation system. Considerable effort has been 
expended in the development of a variety of low- thrust guidance schemes 
(Refs. 1-18) which range from the extremely simple to the mathematically 
elegant. 

However, none of these methods offers a realistic solution to the prob- 
lem of operational guidance when actual mission constraints must be con- 
sidered. In general, all of them generate control programs or control 
program corrections which are unacceptable from an operations viewpoint 
because they require control accelerations that may be beyond the capability 
of the spacecraft (the problem of controllability), or that are continuously 
varying in a way that imposes excessive hardware requirements for imple- 
mentation. In addition, some guidance laws delay corrections until an 
optimum time (with respect to some preselected performance criterion such 
as minimum fuel); and although this policy at first may appear to be desir- 
able, it can lead to controllability problems if additional errors or distur- 
bances occur between the guidance update time and the selected correction 
time . 

The purpose of this memorandum is to present a practical algorithm 
which may be used for both analysis and operational guidance. The develop- 
ment of the steering law is based on the use of linear perturbation theory to 
correct a nominal trajectory in a way that minimizes violations of desired 
terminal mission constraints. The corrections are required to occur prior 
to a prespecified time, thus causing an immediate attempt to improve the 
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trajectory before additional errors occur and possibly give rise to problems 
of controllability. In addition, the implementation deficiencies of earlier 
schemes are avoided by restricting the correction policy to control variable 
changes which are constant for specified time intervals and are of limited 
magnitude. Section II contains a complete derivation of the deterministic 
(i. e. , operational) form of the guidance equations, and Section III presents 
the statistical formulation of those same equations as required for a statis- 
tical guidance analysis. 

In order to evaluate the effectiveness of the basic algorithm, it was 
employed in a preliminary analysis of the approach phase of a rendezvous 
with the comet Encke. Section IV describes that analysis and discusses the 
performance of the guidance scheme, 

II. THE LOW- THRUST GUIDANCE ALGORITHM 

In general, a nominal mission trajectory is designed to satisfy termi- 
nal constraints 


where 



i}j = m vector set of constraint equations 
X (t) = nominal trajectory state, a 6 vector 
t£ = nominal mission final time 

For guidance purposes, it is assumed that Eq. (1) has the form 


where 




0 


y(t) = state of the target set, a 6 vector 


( 1 ) 


(2) 


Since the actual mission trajectory generally deviates from the nominal, 
the terminal constraint equations may not be satisfied at the end of the actual 
trajectory. That is. 
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where 


x(t^) - y(t^) 


0 


(3) 


x(t) = actual trajectory state, a 6 vector 
tj = actual mission final time 
e^ = constraint violations, an m vector 


Equation (3) can be approximated by the linear terms in a Taylor's Series 
expansion about the nominal trajectory: 


4 ^. 


X 






where 

>;< >!< # 

6x^ = x(t^ ) - X (t^ ) = state deviation at t^ 

6tf “ " final time variation 

>;< • >;< 5|< >!s >!< 

x^ = X (t^ ) = rate of change of x (t) at t^ 

• > 1 < • >;< 

yf = y(t^) = rate of change of y(t) at t^ 

4“ = m X 6 matrix of partial derivatives of 4 with respect to 

^ 5j< 

X evaluated at t^ 


(4) 


In order for the actual trajectory to satisfy Eq, (Z), e^ must be zero. 

A natural guidance law, therefore, would be one that generates a control 
program defining a trajectory which nulls out e^. However, because of the 
inherent controllability problems associated with low-thrust flight, it may 
not be possible to eliminate the terminal constraint violations with the avail- 
able control effort. Consequently, a guidance law seeking to null e^ is 
unrealistic. An alternative and more practical approach is to develop a 
strategy which minimizes 


J 


T 

e^Se^ 


(5) 


where 

S = diagonal weighting matrix 


S. . > 0, i = 1 to m 
11 
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A guidance procedure of this type yields the smallest end constraint 
violations attainable with the control available and thus avoids the problem 
of controllability. 


By employing the expanded form (Eq, 4) of e^, minimization of J 
becomes a problem of determining the minimizing final state and time devia- 
tions 6xj and 6t^, Expansion of the actual trajectory about the nominal gives 
6x^ as a function of initial state deviations 6x^ = ^ devia- 

tions from the nominal control program 6u(t) = u(t) - u (t); 


/'■ 


6x^‘ = 4>(t^, t^jax^ + J 4>(t^, T)f^(x, X, u”' )6u(t) dr 

^0 


( 6 ) 


where 


*(t, tj,) 

O.. 

f^(t, X , U ) 

6u(t) 


state transition matrix of the nominal trajectory 

partial derivative of the equations of motion with 
respect to control u evaluated on the nominal trajectory 

control deviations, a 3 vector 


Replacement of 6x^ by Eq. (6) transforms the problem into the minimization 
of J as a function of 6u(t) and 6t^. 

In order to define a practical and useful guidance law, the control 
deviations are required to satisfy a set of auxiliary conditions. Let t be 
a specified time, t^ < t^ , and subdivide the interval (t^, t^) into M equal 
subdivisions (t^, i = 0 to M - 1. Then for 6u(t) to be admissible, 

it must satisfy 

(1) 6u(t) = 6u^ = constant, te(t^, 

(2) 6u(t) = 0, t > t^ 

(3) -a < 6u^ < a, a - specified 3 vector 

Substitution of Eqs. (4) and (6) and the admissible control program 
into Eq, (5) reduces the problem to the minimization of the function 


[ > .X • • V TT r/. • 5l< • >1< -I 

I + r6u + (x^ - y^ )6t^ A I + F6u + (x^ - y^ 
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subject to 


where 


-a < 6u < ct 


A = 

U st^ 

1' n' 

(tf) - 

•j> "1 r 

y(tf )J Svjj^^x 

•X* ^ T 

*if' •'P 1 

(tf) - y(t£)J 

i = «(t* 

to)6xo 





1 '' 1 


^ = ['’o i 

^li" 

• I^M-lJ* 

6 X 3M matrix 


/"‘i+i 

* 

* *v , 

f. = 

1 

L 

1 

^(tf , T )f^(T, 

X , u )dr, 1 = 

^ T Ft 

6u' = 

' ;; T 
1 ou , 
1 t 

' 1 T 

1 1 M- 

3M vector 

Q'1' “ [a'l^ J 


• • 1 YV^I 

,Q J , 

3M vector 


- 1 


Performing the minimization of J with respect to St^ yields 

-1 


f . • >!< T • * • *vl 

6tf = -[(x^ - y^) A(x^ ■ yf)J • 

(Xf - Yf 


) aI + (x^ - y^) AP6u 


( 8 ) 


(9) 


By introduction of Eq. (9) into E.q, (7), J reduces to simply a function of 6u: 

iT, 


= j^l + r6u A 1^6 + r6u 


( 10 ) 


where 


A = 


• s!< •;{< *:{:T 

(x^ - yf)(x^ - Yf) A 


A ' 1 m M. 

(x^ - yf) A(x^ - y^) 


T r 


I - 


• ^ 

(Xf - y£)(x^ - y^) A 

a kI# h A m vt> 

?j> • »-p • *-p 

(X£ - y^) A(x^ - y^) 


A qualitative discussion of the significance of the reduced form (Eq. 10) for 
J is presented in Appendix A. 
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Since J is now a positive semidefinite quadratic function of 6u, it may 
not have a unique or easily computed minimum. These difficulties are 
removed by augmenting J with an additional term to form a new function J: 

J = J + 6u'^W&u (11) 

The matrix W is a positive semidefinite weighting matrix defined by the 
following procedure: 

(1) Find the orthogonal transformation which diagonalizes the 

T-“ 

symmetric matrix F AF, 

Q(f'^AF)q'^ = D 

(f'^AF) = q'^DQ (12) 

where 

D = diagonal matrix, Q = orthogonal matrix 

(2) Define a diagonal matrix B, 

B = diagonal (b^ (13) 

where 

b. = 0 for D. , 0 

1 11 

b. > 0 for D. . = 0 
1 11 

(3) F orm W, 


W = Q^BQ (14) 

T 

Because of the term 6u W6u, J is a positive defiaite quadratic function with 
a unique minimum. Moreover, as shown in Appendix B, the unconstrained 
minimum of J is also an unconstrained minimum of J, and by a suitable 
choice of the b^, the constrained minimum of J closely approaches a con- 
strained minimum of J, 
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Straightforward minimization of J subject to the constraints of Eq. (8) 
yields the feedback control law 

bn = - (r'^AE+ W + A)"^r'^A$(tJ, tg)6xQ (15) 

where 

A = diagonal (\^), > 0 are the Lagrange multipliers which enforce 

the constraints on 6u 

Since the matrix (r^Ar+ W) = q'^{D + B)Q is nonsingular, the inverse 
required in Eq. (15) always exists. Equation (15) can be substituted back 
into Eq. (9) to give 5t^ in feedback form: 






I - r 


T~ 

r AE + 


w+a) 


$(L 




( 16 ) 


The controls 5t^ and ^u are coupled only through the matrix A; consequently, 
in the guidance calculations 5u can be found independently of 6t^, the multi- 
pliers can be determined, and then 6t^ can be computed. 

The basic variable-time-of-arrival guidance algorithm described in 
this section can also be used for fixed-time- of-arrival guidance by setting 
(x* - yj) = 0, A = A, and omitting 6t^ computations. Equation (15) again 
defines the control law 6u, which for this fixed-time case generates a tra- 
jectory minimizing terminal constraint violations e^ at the nominal final 
time t^ , 
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in. STATISTICAL GUIDANCE ANALYSIS 


Section II presented a set of deterministic guidance equations defining 
control corrections based on a given state deviation. In mission planning 
and performance studies, however, statistical analysis techniques are nor- 
mally employed, and only a statistical description of the state deviations — 
the state covariance — is available. Consequently, for analysis purposes, 
it is necessary to have a statistical formulation of the guidance laws which 
utilizes only the state covariance. In a straightforward manner the following 
expressions can be derived for the covariance of the final time and control 
deviations and for the final state errors remaining after the application of 
the guidance strategy: 


r T 

'^‘fj = V = [<==£ - ’'f > *<*£ - 5'f>J <*f - yf> 


[i - r( 


r-^Ar+ w + A 


[i - r(r^Ar+ w +a) e’^aj a(x^ - y^ 


) e'^a]z • 

T 

.] 


(17) 


e[6u6u'^] = u = (r'^AE+ w+a) f'^az A r(r'^Ar + w+a) 


(18) 


E 



X T~ T~1 

= = L 

I - F(F AF+ W +a) F aJ 


I - r(r'^Ar + w +a) f'^a 


where 




- _ r T i 

^0 ~ ^ ^^(tQ)Sx (tp) = state estimate covariance at time t 


(19) 
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The covariances (17 - 19) represent ensemble averages of time, 
control, and final state deviations based on a guidance policy which employs 
the feedback control gain matrix 


G = ~(r^Ar+ w+a) r'^A (2o) 

In the deterministic case, the matrix. A, and therefore G, are selected to 
enforce constraints on the control variables. In a statistical analysis, how- 
ever, only the covariance of the controls is computed, and a question arises 
as to the choice of A. Setting A to zero gives the unconstrained control gain 
matrix and represents an unrealistic choice since unrestricted control is 
unavailable on actual missions. A more reasonable approach is to define A 
such that the standard deviations of the control variables satisfy constraints 
in a manner similar to the deterministic controls, i.e. , 


0 < (Tc £ O' 
6u 


( 21 ) 


where 

or,. = 3M vector of control standard deviations 
ou 

~Q = 3M vector of limiting values 

The use of constraint (21) in selecting A permits practical statistical 
guidance analysis to be performed with a constrained control law. The pro- 
cedure is, therefore, an improvement over previous methods which 
employed the weighting matrix technique to reduce control deviations to 
acceptable levels. 


IV. GUIDANCE ALGORITHM EVALUATION - 
ENCKE RENDEZVOUS GUIDANCE 


In order to evaluate its capabilities, the guidance algorithm was used 
for a statistical analysis of the approach phase of a 1980 rendezvous mission 
with the comet Encke. Orbit determination was accomplished by processing 
daily on-board optical and range observations with a Kalman Filter (Ref. 19). 
The standard deviations of the angle pointing error, target center finding 
error, and range measurement error were taken as 100 arc seconds, 10 km. 
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and 1 km, respectively. The vehicle was assumed to be experiencing a 
random acceleration which was modeled as a first-order Gauss-Markov 
process. The acceleration vector components were assumed uncorrelated 
and spherically distributed with a standard deviation of 1. 8% of the nominal 
thrust-acceleration level and a correlation time of 5 days. At the start of 
the approach phase, the standard deviations of the state errors were chosen 
to represent the relative comet- spacecraft state uncertainty due to the comet 
ephemeris uncertainty. The values selected were 30, 000 km in position 
and 11. 57 m/s (1000 km/day) in velocity. 

The guidance algorithm was employed in four different modes of 
operation; 

(1) Every day a single control correction 5u^ was computed and 
applied for a period of 1 day. The final time was held fixed. 

(2) E very 2 days, two control corrections, Su^ and 6u^, were com- 
puted and applied for successive periods of 1 day each. The 
final time was held fixed. 

(3) Same as (1), but with the final time variable. 

(4) Same as (2), but with the final time variable. 

The control variables for the analysis were the components of the thrust- 
acceleration vector, and the limiting values on their standard deviations 
were set at 10% of the nominal thrust-acceleration level. 

The guidance scheme performance for modes 1 and 2 is shown by 

Figs. 1 and 2, which give the final rms position and velocity errors as a 

function of the weighting matrix S and the time prior to nominal rendezvous 

at which guidance was initiated. For this rendezvous mission, S is a 

6x6 diagonal matrix with the first 3 diagonal elements set to S (the posi- 

P 

tion weight) and the second 3 set to (the velocity weight). Figs. 3 and 4 
give the final rms errors for modes 3 and 4 as a function of S for a 30-day 
guidance initiation time; the corresponding mode 1 and 2 results are repeated 
for comparison purposes. Figs. 5 and 6 show the mode 2 and 4 perform- 
ance when the constraint limit on one of the controls is reduced to 5% of the 
nominal thrust-acceleration; the corresponding results for the original con- 
trol constraints are also included for comparison. 


10 
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The effectiveness of the guidance scheme can be judged simply by 
determining whether the final state errors can be reduced below some 
acceptable limiting values. For the Encke rendezvous, those limits were 
set at 1000 km and 4 m/s. It is clear from Figs. 1 and Z that there exist 
some values of S for which successful rendezvous with modes 1 and 2 is 
possible, provided guidance is initiated 40 days and 35 days, respectively, 
prior to nominal encounter. Moreover, the trends indicated by Figs. 3 and 4 
suggest that for a variable-final-time mode, those initiation times may be 
reduced almost 5 days with little loss in performance. Figs. 5 and 6, on 
the other hand, show that tighter control constraints degrade guidance effec- 
tiveness, and consequently earlier initiation will be required to achieve ren- 
dezvous. Again, permitting final time to vary improves performance. 

General conclusions from the guidance scheme evaluation include the 
following: 

(1) Because trajectory correction with a low-thrust vehicle requires 
a significant time period, early guidance initiation is required 
in order to obtain small final errors. 

(2) The weighting matrix S provides a trade-off in final position and 
velocity errors; the effect of the trade-off is more pronounced 
for the later initiation times, which have the greater overall 
final errors. 

(3) The two- correction policy is better than the single-correction 
one for rendezvous, since 6 rather than 3 control variables are 
being used to control the 6 final state deviations. But for a flyby 
(S = 1, S = 0), two corrections offer no improvement over 

' p V 

one. 

(4) The variable-final- time mode of operation significantly improves 
performance by providing an additional control variable. 
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V. CONCLUDING REMARKS 


This memorandum described and evaluated a low-thrust guidance 
algorithm that has been designed to allow for actual mission constraints. 

It is a linear perturbation scheme with the following primary features: 

(1) Violations in the terminal mission constraints are minimized. 

(2) Control corrections are required to occur early in the trajectory. 

(3) Control corrections are constant over specified time intervals. 

(4) Control corrections are explicitly bounded. 

By having these features, the algorithm avoids the controllability and imple- 
mentation difficulties associated with other guidance laws. 


12 
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Fig. 2. RMS final state errors vs flight time and 
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APPENDIX A 

THE REDUCED J FUNCTION AND THE FINAL TIME VARIATION 
The function J (Eq. 10) can be written in the form 


J = (6x^)^A(6x^) 


where 


A = 


T 

E AE 


E = I 


-[<*£ 


• *.T. 

Yf) A(x^ - y^) (x. 


si^ 5k 


y£)(xj - Yf) A 


Let the time derivative of the terminal constraints, evaluated at the end of 
the nominal trajectory, be denoted by the vector q, 


q = + 


X 


[ Kt» vO I 

X (t£) - y(t£)J(A£ - y^) 


and let the vector p be defined as 


[>!< ❖ ❖ n ❖ 

X (t^) - y(t^ )jE6x^ 


In a straightforward manner, it may be shown that when S is the identity, 
p is orthogonal to q: 


i .*.Tj^T 

q P = (Xf - y.) V 


f ^ f ' 'X 


x'‘'(tj') - y{t;')j4^^|x^'(t;) - y(t;)jE6x; 


» 'T* « •'I'* I 'V* ^ 

- y^) AE5x^ = 0-6x^ = 0 


The vector p represents the constraint violations at t^ which are orthogonal 
to the nominal time derivative of the constraints. For a suitable choice of 
unitsfor e^, the matrixSmay always be selected as the identity; consequently, 
it is clear that minimization of the reduced form of J as a function of 6u 
simply yields the smallest possible p, i.e. , J = p'^p. The final time varia- 
tion 6t^ is then used to null out remaining terminal errors along the direc- 
tion q. 
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APPENDIX B 

THE WEIGHTING MATRIX W 


The primary functions of the weighting matrix defined by Eqs. (12-14) 

are 


(1) To guarantee a unique minimizing control 6u, 

(2) To permit the minimizing control to be computed in a straight- 
forward and systematic manner. 

(3) To guarantee that the control 6u can always be expressed in feed- 
back form, i.e., Eq. (15). 

For the particular choice of W, it can be shown that in the unconstrained 
case, the control which minimizes the function J (Eq. 11) also minimizes 
the function J (Eq. 10). In order to verify this property, the following 
lemma is required. 

LEMMA: Let the m x m orthogonal matrix Q be the diagonalizing transfor- 

mation for the matrix (P AP), where A is an n x n positive semidefinite sym- 

T~ T 

metric matrix and P is an arbitrary n x m matrix, i, e. , QP APQ = D. 

T 

The matrix D is an m x m diagonal. Let q = QP A§ be an m vector where 
% is an arbitrary n vector. Then if element D_ of matrix D is zero, com- 
ponent of q is also zero. 

Proof: Let the orthogonal matrix V diagonalize A, then VAV^ = A, where 

A is diagonal. It follows that (QP^v'^)A(VPQ^) = D. Let p.. be an element 
of matrix Qp'^v'^, then 


n 

y p^.A.. = D.. 

Z. H JJ 11 

j = l 

for each i. 

Since A is positive semidefinite. A.. > 0. Consequently, if D.. = 0, 

2— , _ JJ ^ ^ 11 

then p..A.. = 0 for all j, or p..A.. = 0 for all j. The vector q may be 
JJ T T — 

written as q = (QP V )A{Y%). The i'th component then is 


^ _ 

q. > p..A..r. 

1 ZL IJ JJ J 

j = l 
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where r. is a component of the vector r = V^, If D.. = 0, p..A.. - 0 for 
J I 

all j; consequently, = 0 as required. 

With the use of the above lemma, it is easy to establish the following 
property of the unconstrained control minimizing J, 

Property: The control 6u which minimizes the function 

J = [^+ P6u]'^a[^+ P6u] + 6u'^W6u 

also minimizes 

J = [|+ r6u]'^A[|+ P6u] 

Proof: Let the vector w be defined as w = Q6u + (D + B) then 

^ T n X n 

6u = Q'*'[w - (D + B) Substitution of this expression into the function J 

yields 

J = ^"^Al + w'^Dw + 2w^[l - D(D + B)"^]q 
-q'^(D + B)"^[2(D + B) -D](D + B)'% 

If D.. = 0, it follows by definition that B.. >0, and by lemma that q. = 0; 
also, if > 0, B^^ = 0. Consequently, J may be reduced to the form 

J = w^Dw + C'^AI - q'^(D + B)'^D(D + B)'% 

A minimum of J clearly occurs at w = 0; that is, Q6u + (D + B) ^q = 0, 

^ T - 1 

which implies 6u = -Q (D + B) q. But this control is the one obtained by 
minimizing J; therefore, the control which minimizes J also minimizes J. 

In the constrained case, 6u may not possess the above property. How- 
ever, by selecting the elements of B sufficiently small when forming W, the 
minimum of J can be approached to any desired accuracy. Let 6u mini- 
mize J subject to the constraint - S’ < 6u < 3, then it is necessary that 

(P^AP + W)6u + P'^A^ = -A6u (B-1) 
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where A is the diagonal matrix of non-negative Lagrange multipliers. For 
B sufficiently small, it is clear there exists a non-negative diagonal matrix 
A such that 


A6u + W6u = A5u + c (B-2) 

where e is a small residual error vector. Therefore, by substitution of 
Eq. (B-2) into Eq. (B-1), 


rp rp 

(F AF)6u+ F A4 «=> -A6u 

but this is simply the condition for the minimization of J. Consequently, 

'6u is an approximation to the constrained minimizing control for J. For 

computational purposes, the elements of B should be small with respect to 

those of D, but should be large enough to avoid numerical difficulties in the 

T~ - 1 

computation of the inverse matrix (F AF + W) 
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